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ABSTRACT 

We investigate the radio emission of ~ 185,000 quiescent (optically unclassifiable) galaxies selected 
from the Sloan Digital Sky Survey (SDSS). By median- stacking FIRST cutouts centered on the 
optically-selected sources, we are able to reach flux densities down to the 10s of /xJy. The quies- 
cent galaxy sample is composed of two subgroups inhabiting vastly different regimes: those targeted 
for the SDSS MAIN Galaxy Sample (~55%), and those targeted for the Luminous Red Galaxy (LRG) 
sample (~45%). To investigate the star- format ion of these quiescent galaxies, we calibrate a radio-SFR 
conversion using a third sample of star- forming galaxies. We confirm a tight power-law dependence 
for the star-forming (SF) sample, where Liaghz ~ (SFR) 1 - 37 . Comparing this SFR-indicator with 
indicators in the optical and UV, we derive conflicting SFR estimates for the MAIN sample quiescent 
galaxies. These radio-derived SFRs intersect those calculated using the 4000-A break (D4000) around 
an SFR of 1 M© yr _1 and agree to within a factor of 3 over the range of SFRs. However, we find that 
the radio-derived SFRs are too high relative to the SFRs estimated for similar populations of galaxies 
using analysis of UV emission, implying either contamination of the radio by Active Galactic Nuclei 
(AGN) or incomplete dust modeling. If AGN activity is dominant in these galaxies, then a relation 
between AGN radio luminosity and galaxy mass is required to explain the observed trends. For the 
LRGs, on the other hand, we find the radio luminosity to be relatively high (compared to the SF 
galaxies) and independent of SFR as derived from D4000, indicating an AGN component dominates 
their radio emission. AGN-based radio emission often implies the existence of radio jets, providing 
evidence of a mechanism for low-level feedback in these quiescent LRGs. 
Subject headings: 



1. INTRODUCTION 

For normal galaxies - those not deriving their radio lu- 
minosity from a nuclear "monster" — the radio em ission 
is a probe of recent star- formation (| CondonI 1 1992T ). In 
particular, the emission below 30 GHz is thought to be 
largely dominated by synchrotron radiation from ultra- 
relativistic electrons accelerated in supernova remnants. 
As most of the remnants responsible for this non-thermal 
emission result from massive, short-lived stars, the result- 
ing radio emission tracks recent star-formation. Further 
evidence in favor of this claim comes from the existence 
of a tight radio/far- infrared (FIR) correlation, which has 
been found to h old over several orders of magnitude in ra- 
dio luminosity (lHarwit fc Pacinilll975l : ICram et alj|199§ 



lYun et alJl200H 

The complex astrophysics connecting the integrated 
radio emission to supernova explosions has led to a his- 
torical preference for other methods. Normal galaxies, 
with quite low radio flux densities, are also usually be- 
low the sensitivity of typical wide-area radio surveys. Al- 
ternative indicators of star-formation rate (SFR) include 



continuum e mission in both the near-ult r aviolet and the 
far-in frared (|Donas fc Deharvend Il984t lHarper fc Low! 
1 19T3T ) . Another common diagnostic, referred to later in 
this paper, is the use of integ rated emission- line fluxes 
(|Cohenill976t lKennicuttl[i983li . (For additional indica- 
tors, see the excellent review bv iKennic utt (1998)).) 

However, star-formation is not the only process that 
produces radio emission; Active Galactic Nuclei (AGN) 
can also be strong emitters in the radio. This ambiguity 
has led to much debate, particularly for the less well- 
studied sub-mjy population. On the one hand, there is 
a widely held belief that the population of sub-m jy radio 
sourc e s is dominated by star - forming galaxies (ICondonl 

iKron. Koo. fc Windhorst! 
For e x ample , using an 
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infrared-selected sample. iBovle et all (|2007|) found the 
radio/FIR correlation persisted down to /zJy flux den- 
sities, and they conclude that the radio emission of 
these sources is dominated by star-formation (with red- 
shift still unknown). However, the hypothesis that there 
are a significant number of AGN at these low flux 
densities has recently re ceived some attention as well. 
Uarvis fc Rawlingj (|2004f ). in a theoretical paper, sug- 
gested that the flattening in the observed counts of the 
normalized local luminosity function below S ~ 1 mJy 
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is due to a "radio-quiet " AGN popula t ion. I n a follow- 
up observational study, ISimpson et al.1 (|2006l ) concluded 
that radio-quiet AGN comprise a significant population 
below ~ 300 /xJy, and that they may, in fact, dominate. 

Studying the general properties of the low-luminosity 
radio population, however, has proven to be difficult; the 
perennial trade-off between coverage area and depth dis- 
courages surveys that are both deep and wide. In this 
paper, we take advantage of the overlapping coverage of 
two surveys, FIRST (Faint Images of the Radio Sky at 
Twenty Centimeters) and the SDSS (Sloan Digital Sky 
Survey), to study the faint radio emission of a large sam- 
ple of optically- normal ("quiescent") galaxies. By select- 
ing galaxies with no strong optical emission lines from 
star-formation or AGN activity, we construct a sample 
that is likely to be inherently radio-quiet in nature. Any 
radio emission we do find, then, must result from either 
extremely low-level SF and/or AGN activity, or activity 
that is obscured in the optical. By observing how this 
emission behaves as a function of other SFR indicators, 
we will attempt to discern the origin (SF or AGN) of 
the radio flux, and, in the absence of AGN activity, the 
radio emission will give us a handle on the amount of 
st ar-formation occur ring in these galaxies. 

Sali m et al.l (|2007f ) found that galaxies with no Ha 
detectable in their SDSS spectra were mostly normal 
early-type galaxies with predominantly red optical col- 
ors, though some have relatively blue UV to optical 
colors (see also: Rich et al. 2005). They calculated 
that some of these galaxies had UV-derived SFRs as 
high as 1 M yr -1 , and a significant portion of those 
showed structure, like traces of faint spiral arms, or dis- 
turbed light profiles. They identified an evolutionary 
sequence from star-formation to quiescence via nuclear 
activity. Using a sample of nearby, early-type galaxies 
in the MOSES project (Morphologically Selected Ellipti- 
cals in SDSS) . TSchawinski et al.l (l2007f) reached a sim- 
ilar conclusion. lYan et alj (|2006l ) looked at optically- 
quiescent galaxies with red rest-frame colors and found 
that they were virtually indistinguishable in the color- 
magnitude-concentration space from galaxies classified 
as low- ionization nuclear emission- line regions (LINERs). 
Still, very little has been done on the radio properties of 
optically-quiescent galaxies, and we intend to use this 
relatively unexplored wavelength regime to add to our 
knowledge of such galaxies. 

For the analysis of such faint radio sour ces, we have 
utilize d the stacking technique presented in I White et al.l 
(|2007f ). hereafter referred to as Paper I. This technique 
significantly reduces the noise in the combined images, 
allowing the investigation of median radio properties for 
sources significantly below the FIRST detection thresh- 
old without requiring countless more hours of telescope 
time. In addition, stacking all of the fields guarantees 
that the radio data will not be restricted by survey lim- 
its. In de Vries et al. (2007; hereafter, Paper II) we ap- 
plied this method to the analysis of low- luminosity AGN, 
concluding that radio emission from star-formation can 
dominate over that due to the AGN emission in compos- 
ite systems. 

We begin by describing the radio data and the details 
of the sample selection ($2j. We then go into the details 
of the stacking method ($3} • In S2) we present our analysis 
and a detailed discussion of the results; §b\ contains our 



concluding remarks. Throughout this paper we use the 
standard FRW cosmology of Ho = 70 km s -1 Mpc -1 , fl\ 
= 0.7, and Q M = 0.3. 

2. SAMPLE CONSTRUCTION 
2.1. The Radio Data 

The 1.4 GHz radio d ata we present are taken from th e 
VLA FIRST Survey (|Becker. White, fc Helfandl fl995T) . 
which includes 9,030 deg 2 sky coverage of the North 
Galactic Cap and Equatorial Strip, making it ideal for 
comparison with targets from the Sloan Digital Sky Sur- 
vey (SDSS). With only - 90 sources deg -2 at the 5cr 
detection limit, 99.98% of the 5 billion beam areas are 
'blank' (Paper I). The data, obtained from 180 sec VLA 
snapshots, have a typical rms of 0.15 mJy beam -1 . By 
stacking 3000 such snapshots, we can reduce the noise in 
the resultant image to approximately 3 /iJy beam -1 and 
achieve the sensitivity equivalent of 150 observing hours. 
To convert to rest- frame 1.4 GHz luminosity, we perform 
a k-correction assuming an average radio spectral index 
a = -0.5. This index provides a compromise between 
steep-spectrum radio-lobe emission and flatter spectrum 
emission from the AGN cores or star-forming regions, 
and it is the same index adopted throughout Paper II. 

2.2. The Quiescent Galaxy Sample 

The sample of quiescent galaxies w as drawn from the 
SFR Catalog of the MPA/ JHU VACs (jBrinchmann et all 
l2004h . The SFR catalog is a super-s et of the galaxies 
discussed in iBrinchmann et all (|2004h and contains all 
galaxies in the SDSS DR4 spectroscopic survey. DR4 of 
the SDSS has 6851 deg 2 of i maging and 4681 deg 2 of 
spectr oscopic coverage. (See lAdelman-McCarthv et al.l 
(2006) for more information on DR4). Our final sam- 
ple of quiescent galaxies is complete to a stellar mass 
lower limit of 10 8 M . To accurately classify the ob- 
jects, Brinchmann et al. reanalyzed the 1-D SDSS spec- 
tra using their own o ptimized pipeline as outlined in 
iTremonti et al.l (|2004h . A galaxy was classified as ei- 
ther star-forming (SF), low S/N SF, AGN, low S/N 
AGN, Composite (C), or Unclassifiable based on the line 
ratios of the four emission lines described in Baldwin, 
Phillips, & Terlevich (1981, hereafter BPT). Composite 
galaxies have line ratios intermediate between SF and 
AGN. For the SF and low S/N SF classes, the fiber 
SFR was obtained from detailed modeling of the emission 
lines, while for the AGN, C, and Unclassifiable galax- 
ies, the fiber SFR was estimated by calibrating the re- 
lation between the specific SFR (calculated for the star- 
forming gakxies}_and the m easured D4000 value as de- 
fined in lBalogh et~"aTl ([1998). We will refer to the latter 
as "D4000" SFRs to differentiate them from SFRs ob- 
tained via emission lines. For both types, total SFRs 
were calculated using an empirically based aperture cor- 
rection scheme based on the likelihood distribution of 
the specific star formatio n rate for a given set of colors 
(|Brinchmann et al.ll2004h . 

To single out only "quiescent" galaxies, we kept only 
those galaxies catalogued as Unclassifiable. The majority 
of thi s class are galaxies with no or very weak emission 
lines (jBrinchmann et al.ll2004l ). As the SDSS uses fiber- 
fed spectrographs, this is a 'nuclear' classification. The 
resultant sample contains ~ 185,000 quiescent galaxies. 
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Fig. 1. — A median image resulting from stacking 3000 quiescent 
LRGs. The image is 0.5' X 0.5', (larger than normal, so the reader 
may observe the obvious sidelobes), and each pixel is 1.8" on a 
side. The cutouts comprising this image corres pon d to the galax- 
ies making up the right-most circle in Figure 1101 prior to being 
converted to radio luminosity. The image has a peak flux density 
of 90 /^Jy, with an rms of 6 /ijy beam -1 . This image is virtually 
indistinguishable from any of the other high S/N images generated 
for this study, independent of the type of galaxy. 

We obtain additional parameters of interest for each 
galaxy from the NYU Value-A dded Galaxy Catalo g 
(NYU-VAGC). As presented by iBlanton et aD ((20051) . 
the NYU-VAGC is a collection of galaxy catalogs derived 
from the SDSS. The latest version of the catalog accom- 
panies the SDSS DR4. We also make use of the derived 
ugriz absolute Petrosian magnitud es available in the as - 
sociated catalog of K-corrections (|Blanton et al.l [2007^ . 
As these data were originally computed using a different 
value for the Hubble constant (Ho = 100), we adjusted 
the derived absolute magnitudes to agree with the stan- 
dard cosmology used elsewhere in this paper. 

2.3. The Star-forming Galaxy and AGN Samples 

The star-forming sample, first introduced in Paper II, 
consists of all galaxies classified as star-forming (SF) 
in the DR2 SFR catalog. This sample has S/N > 
3 in all four emission lines of the BPT diagram and 
falls below the conservative threshold for AGN rejec- 
tio n. Star-formation rates quoted here were derived 
by iBrinchmann et aTl ((2004) via careful modeling of the 
emission lines, and they take into account dust attenua- 
tion and fiber-c orrections. No te that these SFRs assume 
a Kroupa IMF (|Kroupall200lD . For absolute magnitudes 
and derived parameters , we consulted the Stellar Mass 
Catalog as presented bv lKauffmann et ail (|2003bl) . 

The AGN sample was taken directly from Pap er II and 
consists of objects from the iKauffmann et all |2003a) 
AGN compilation. It should be noted that they use a 
more liberal AGN cut than IBrinchmann et aT] (|2004f) . 
and this sample therefore consists of objects classified as 
AGN as well as those classified as Composite in Binch- 
mann et al. This is a simple naming-convention issue and 
does not affect our results. It is important to remember, 
however, that objects in the AGN sample are affected by 
star-formation to different degrees, with those closest to 
the cut having emission-line ratios more commonly asso- 
ciated with HII regions. We use the D4000 SFRs (as we 
do wi th the quiescent galaxies) from Brinchmann ct al. 
(200|). 



3. STACKING FIRST IMAGES 

Using the sample selection criteria described above, we 
compiled our optically-selected quiescent galaxy sample. 
We then obtained 1.4 GHz data centered on each target 
position from the FIRST survey database. These square 
"cutouts" of the radio sky were 0.2' on a side and had the 
typical FIRST rms of S rm s = 0.15 mJy beam -1 ; these 
cutouts constitute our 3-D "image block" . The size of the 
images was chosen pragmatically based on limitations of 
IDL memory. Since all the analysis is based on the value 
of the central pixel, the image size is not important. In 
any case, all stacked images have a FWHM < 8" and 
hence easily fit within the image format. These blocks 
can be sorted by a variety of parameters (e.g., redshift, 
color, SFR, etc.) and stacked to yield median radio flux 
density/luminosity as a function of the parameter of in- 
terest, with a bin size (N bj / bin) determined by the 

radio survey depth required (S m i n ~ S rms / y/~N). 

Bins typically contained either 1500 or 3000 images 
each. The median of the parameter of interest was de- 
termined for each bin. To obtain the associated radio 
luminosity, the cutouts were processed through an IDL 
procedure which converted each individual image to units 
of radio luminosity using the target galaxy's measured 
redshift (Paper I) and then output the median image 
and ±lcr upper/lower bounds for each pixel. The basic 
algorithm for determining the error bars involved com- 
puting the probability that the true median lay between 
elements (k,k+l) in the sorted i mage stack using bino- 
mial probabilities as described in lGott et alj (|2001h . As 
these sources are all quite faint in the radio, we will sim- 
ply quote the radio luminosity derived from the central 
pixel in each co-added image as done in Paper I. 

The median was used as opposed to the mean because 
it is more statistically robust in the presence of outliers 
and/or contamination. For example, S alim et alj (|2007f ) 
note that ~ 2% of the unclassifiable galaxies have line 
emission, some quite significant. The median will be lit- 
tle changed by a 2% contaminant. The mean, on the 
other hand, is sensitive not only to outliers, but also to 
the threshold values imposed to try to weed out those 
outliers (Paper I). When we used means, we found the 
radio luminosities all increased by 1-2 orders of magni- 
tude. That being said, the interpretation of the median 
is slightly more complicated. The median shifts toward 
the 'local' mean in the case of low S/N data. (See Paper 
I for a more detailed discussion). 

Once the radio data were stacked, we corrected for the 
so-called 'snapshot bias'. This phenomenon is somewhat 
analogous to the 'CLEAN bias' reported for the FIRST 
and NVSS surveys, where the result of the non-linear 
CLEAN algorithm is that a small amount of flux from 
t he above-threshold sources is redistributed over the f ield 
(|Becker. White, fc Helfanc3ll995UCondon et al.lll998f) . A 
similar bias for sub-threshold sources (that have there- 
fore never been CLEANed) was discovered and calibrated 
in Paper I. For sources fainter than 0.75 mJy, only 71% 
of the true mean flux density is recovered per bin, while 
sources above 0.75 mJy are well-described by the ad- 
dition of a constant 0.25 mJy CLEAN bias. Since the 
majority of the target galaxies in this study are sub- 
threshold, we accounted for this snapshot bias via the 
constant multiplicative factor of 1.40. 
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Density vs. Apparent Magnitude 




i i 



Fig. 2. — A plot of radio flux density versus apparent magni- 
tude for the entire quiescent galaxy sample. The median-stacking 
method allows us to probe flux density levels down to the 10s of 
/ijy. The clear selection bias toward stronger radio sources is evi- 
dent below m z ~ 17; this is a result of the SDSS target selection 
algorithm (see text). 



Radio Luminosity vs. Optically-derived SFR 
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FlG . 4. — Radio Luminosity (1.4 GHz) vs. optically-derived SFR 
from iBrinchmann ct al. (2004). This plot shows the star- forming 
galaxies and AGN originally presented in Paper II, as well as the 
quiescent MAIN sample galaxies and LRGs overplotted. Note that 
all of these objects are still at least 3 orders of magnitude fainter 
than the typical radio-loud object. 



Flux Density vs. Apparent Magnitude 
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Fig. 3. — A plot of radio flux density versus apparent magnitude 
for the quiescent galaxy sample, this time with galaxies targeted 
as part of the MAIN Galaxy Sample plotted separately from those 
targeted as LRGs. The LRGs have much higher median flux den- 
sities. 

An example of a median-stacked image is shown in 
Figure [TJ The image is 0.5' on a side and resulted from 
stacking 3000 quiescent LRGs. The cutouts that went 
into this image were the flux density maps of the galax- 
ies that were later converted to units of luminosity and 
stacked to give the right-most quiescent LRG bin (circle) 
in Figure [TUJ The image peaks at 90 /iJy and has an rms 
of 6 /zJy beam -1 . 

4. RESULTS 

4.1. Quiescent Galaxy Sample Selection Effects 

In order to explore the radio emission of the quiescent 
galaxy sample, we used the median-stacking technique 
detailed in <JJ] Typical radio flux density levels attain- 
able via this method are shown in the plot of flux density 
vs. apparent magnitude (z-band) in FigureO Each point 
is the result of median-stacking 3000 quiescent galaxies. 
As Figure demonstrates, a pronounced selection effect 
becomes apparent when the data are represented in this 
manner. This effect is the result of the SDSS target se- 
lection algorithm which includes different samples; the 
MAIN Galaxy Sample, the Luminous Red Galaxy (LRG) 



sample, and the Quas ar, Star, or Serendipity samples 
(jStoughton et al J 120 02). Galaxy targets are selected for 
inclusion in the SDSS MAIN galaxy sample if they satisfy 
a num ber of criteria outlined in detail in I Strauss et al.l 
(2002). Approximately 55% of our sample of quiescent 
galaxies were targeted as part of the MAIN Galaxy sam- 
ple. The remaining 45% of our quiescent sample, save a 
small percentage of miscellaneous objects which we dis- 
carded, were targeted as part of the LRG sample. As 
described by lEiscnstci n et ail ((2001) , the LRG sample 
extends to higher z and to fainter magnitudes than the 
MAIN galaxy sample and consists of objects that are 
luminous (L > 3L») and intrinsically very red. 

Objects may also have more than one target flag set 
at a time. Of particular concern to us are those objects 
targeted for both the MAIN and LRG samples. Since the 
LRG group is generally more restrictive, and for reasons 
that will become apparent in our 21 we have chosen to 
keep these objects with the L RGs rather than the M AIN 
galaxy sample. As advised in lEisenstein et aTl (|2001l ). we 
then keep only those objects with the GALAXY_RED 
flag set and redshifts z > 0.2. These further restrictions 
account for the breakdown of the luminosity cut at lower 
redshifts and ensure the purity of the LRG sample. 

If the data are first segregated into MAIN galaxies 
and LRGs, a rather different picture emerges. Figure [3] 
again shows flux density vs. z-band apparent magnitude, 
but this time the sources targeted as part of the MAIN 
Galaxy Sample are plotted as diamonds, while those tar- 
geted as part of the LRG sample are plotted as crosses. 
The MAIN Galaxy Sample show typical radio flux den- 
sities in the 10s of /^Jy, with an increasing number of 
Southern Equatorial Strip sources comprising the low- 
magnitude tail. Of these galaxies, fewer than 0.1% were 
found to be coincident with cataloged FIRST sources. 
The LRGs have significantly higher median radio flux 
density levels (though still far below the FIRST thresh- 
old for detection) and exhibit a trend toward lower flux 
densities with increasing magnitude. Out of the LRG 
sample, 0.9% are found to be within 2" of a source in the 
FIRST catalog, a detection rate nearly ten times that of 
the MAIN Galaxy Sample. 
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4.2. Radio Luminosity as a Proxy for Star- Formation 

As mentioned in the Introduction, radio emission can 
be used as a tracer of star-formation for galaxies with no 
contaminating nuclear activity. We use the star-forming 
galaxy sample, described in Section [21 to obtain the re- 
lationship between star-formation rate and radio lumi- 
nosity. The SF galaxies have measured SFRs based on 
Ha and other optical emission lines, and we sorted the 
sample by this "Ha-derived" SFR and binned/stacked 
the data using the method outlined in Section [3l The 
result, published in Paper II and repeated in Figure 0] 
shows a tight correlation between radio luminosity and 
SFR for the SF galaxies which extends over more than 
two orders of magnitude in SFR. The least-squares fit 
to the data gives log(Li 4 ) = (1.37 ± 0.02) log(SFR) + 
(27.67 ± 0.01). 

There have been a number of studies deriving SFR 
from 1.4 GHz emission, m any of which u se the radio- 
SFR relation presented by ICondonl (11992ft . This rela- 
tion, and those like it (|Yun et alj I2001D assume a lin- 
ear relationship between SFR and radio luminosity. In- 
deed, the tightness of the radio-FIR correlation was first 
interpreted as evidence for a linear relationship. How- 
ever, several authors (|Price k. Durid Il992t [Belli 12003ft 
have shown that the synchrotron component, which gen- 
erates ~ 90% of the emission at 1.4 GHz, depends non- 
linear ly on the SFR. iPrice fc Duricl ([19921 ) found that, 
by decomposing the radio emission in the radio-FIR re- 
lation into thermal bremsstrahlung and non-thermal syn- 
chrotron components, the non-thermal radio luminosity 
scaled as (SFR) 1 - 2 , even while the t h ermal emission was 
directly proportional to SFR. iBelll (|2003l ) argued that 
the radio-FIR correlation is a conspiracy resulting from 
the non-linearity of both SFR indicators. They find that 
Lia goes as (SFR) 13 , a result which is very close to our 
derived power-law dependance of 1.37. 

4.3. Radio Emission from Quiescent Galaxies 

Figure [4] shows our results for the quiescent galaxies 
overplotted on the SF galaxies and AGN of Paper II. The 
first thing to notice is the large difference between the be- 
havior of the MAIN and LRG-targetcd samples. Keeping 
in mind that the LRG bins do have significantly higher 
median redshifts than the MAIN galaxy bins, the sam- 
ples nevertheless exhibit very different trends. The LRG 
sample maintains relatively constant median radio lumi- 
nosity even with optically derived SFRs (from D4000) 
varying from 0.1 to 10 M Q yr -1 , suggesting that the ra- 
dio emission from these optically-normal LRGs is domi- 
nated by AGN activity rather than star- formation. I.e., 
if some of the radio luminosity came from star-formation, 
there would be a dependance on SFR. Indeed, for LRGs 
with SFR < 1 M Q yr" 1 , less than 1% of the radio flux 
should originate from SF activity. It has been pointed 
out to us that perhaps the SFRs in the MPA/JHU cata- 
log were never intended to be appropriate for the LRGs 
in the catalog, seeing as they derived from star-forming 
galaxy colors and reach quite spectacular rates for galax- 
ies believed to be red and dead. This does not change 
our result. Whether the LRG SFRs are in the ballpark 
or complete nonsense, the median radio luminosity of the 
sample is simply too high to be attributable to SF alone. 

One may ask how LRGs actually classified as AGN 



compare. We therefore compiled a sample of galaxies 
classified as AGN or Composite from the MPA/JHU SFR 
catalog that were also targeted as LRGs. We found only 
~ 600 such galaxies, not even enough for a full bin and 
not nearly enough to examine trends in the stacked be- 
havior. The partial bin we were able to make lies basi- 
cally right on top of the quiescent LRGs, with a median 
SFR of ~ 4 M yr" 1 and an L B , of 10 29 9 ergs s" 1 Hz" 1 . 

Figure [4] therefore provides evidence that quiescent 
LRGs host AGN, a result that has some interesting im- 
plications. Most current models of cosmic structure for- 
mation include AGN feedback, and both radiative and 
mechanical heating of the "radio mode" AGN are now 
recognized as important in galaxy formation scenarios 
(for e.g., Croton et al., 2006, and Granato et al., 2004). 
As radio emission from black holes often indicates the 
presence of jets, these data lend observational support 
to the existence of a mechanism for low-level feedback 
in these galaxies. In addition, this result also gives a 
quantitative feel for the feedback subsequently implied. 
The constant median value of L# ~ 4 x 10 29 ergs s _1 
Hz -1 radiated by these quiescent LRGs, while still tech- 
nically low-luminosity, is over two orders of magnitude 
larger than that attributed to the "ground state" of the 
optically-detected low-luminosity AGN in Paper II. As to 
why the LRGs are so much more radio-luminous, perhaps 
this is evidence that more massive AGNs are more radio- 
luminous. But speculation aside, it is true that here we 
are dealing with median statistics and therefore have no 
knowledge of the underlying luminosity distribution (i.e., 
whether it is bimodal, etc.). We also cannot say whether 
the spread in the distribution is due to time variability or 
static differences between galaxies. However, by treating 
this value of radio luminosity as a time-average (rather 
than ensemble-average), we get an estimation of the av- 
erage luminosity of radio jets which may result in AGN 
feedback in LRGs. 

We now turn our attention to the MAIN sample galax- 
ies in Figure[4] These galaxies show a trend in radio lumi- 
nosity with SFR which intersects the SF sample around 
1 M Q yr -1 . The agreement is worse for higher values of 
SFR, where the MAIN sample galaxies deviate system- 
atically from the SF galaxies by a little over a factor of 3. 
There are two competing theories for the origin of the ra- 
dio emission from these galaxies: that it is due to SF, or 
that it is due to AGN activity. The fact that \jr shows a 
dependence on D4000 SFR, and that it disagrees with the 
SF sample by no more than a factor of 3 over the whole 
range of SFRs, suggests that it may arise in part from 
star-formation. Alternatively, one may argue we are see- 
ing the results of a correlation between Lr and mass that 
is due to AGN activity, a valid argument that we elab- 
orate on below. However, since the Lr is near the level 
of the SF galaxies, we think it only fair to entertain the 
notion that it is due to SF. If we assume this to be the 
case, the disagreement between the MAIN sample and 
SF galaxies may be an artifact of the different methods 
(D400 vs. emission lines) used by iBrinchmann et aTl 
(2004) to measure SFRs f o r the different samples. Al- 
though IBrinchmann et aTl (|2004f ) proved the reliability 
of these D4000 values for galaxies in the star-forming re- 
gions of the BPT diagram, that doesn't guarantee that 
the calibration yields sp ot-on values for ga laxies with 
suppressed emission lines. ISalim et"aTI 1)20071 ) argue that 
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Quiescent MAIN Sample: Composite Spectrum for Lowest SFR Bin 



3000 4000 



1/ fAi 




3650 3700 3750 3800 6500 6550 6600 6650 



5000 6000 7000 

Wavelength (A) 



Quiescent LRG Sample: Composite Spectrum for Lowest SFR Bin 




3650 3700 3750 3800 6500 6550 6600 6650 



3000 4000 5000 6000 7000 

Wavelength (A) 



Quiescent MAIN Sample: Composite Spectrum for Highest SFR Bin 
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Fig. 5. — Composite spectra for the MAIN sample quiescent 
galaxies. The top panel (a) shows the composite spectrum for the 
3000 galaxies comprising the bin with the lowest SFR in Figure [J] 
while the bottom panel (b) shows that of the bin with the highest 
SFR. The insets show the spectral features associated with Oil, 
Ha, and NIL 

the D4000 - specific SFR calibration and the color-based 
aperture corrections are only completely accurate for the 
SF galaxies, since both methods are dependent on time- 
scale and other galaxy types have different star-formation 
histories. 

The wary critic might ask whether D4000 is measur- 
ing star-formation at all in this case, or whether AGN 
emission is contaminating the D4000 break as well, lead- 
ing to overestimation of the D4000 SFRs for the MAIN 
sample. We have considerable reason to believe that 
this, at least, is not the case. For one thing, the ab- 
sence of observable line emission from the AGN makes 
it unlikely that the AGN would contribute significantly 
to the continuum near the 4000 A br eak. And further- 
more, e ven for optically-detected AGN.lKauffmann et al.l 
(|2003af) showed that the AGN contribution to the blue 
continuum is ignorable. 

To get a better feel for these objects, we take ad- 
vantage of our large sample size and examine trends in 
median-stacked spectra. We have created median com- 
posite SDSS spectra for both the MAIN galaxy sample 
and the LRGs as a function of SFR bin as shown in 
dJ These spectra, normalized such that the continuum 
value at 4000 A is unity, are shown in Figures [5] and [6l 
respectively. In each case, the top panel shows the spec- 
trum for the lowest SFR bin and the bottom shows that 
of the highest SFR bin. We should note that we also 
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Fig. 6. — Composite spectra for the quiescent LRG sample. The 
top panel (a) shows the composite spectrum for all the galaxies 
comprising the bin with the lowest SFR in Figure [J] as calculated 
from D4000, while the bottom panel (b) shows that of the bin with 
the highest SFR. The insets show the spectral features associated 
with OIL Ha, and NIL 
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Fig. 7. — [Oil] Luminosity vs. SFR for the quiescent MAIN and 
star-forming galaxy samples. Values for [Oil] Luminosity are rough 
es timates using the conversion from equivalent width to luminosity 
in lKenn icutt (1992) with median g'-band luminosity. 

tried binning spectra by SDSS z-band absolute magni- 
tude, and the resulting composite spectra had line emis- 
sion no more obvious than those spectra binned by SFR 
that we show here. 

These high S/N composite spectra clearly show very 
weak emission features. They do, however, show strong 
Balmer line absorption; strong enough to plausibly swal- 
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Fig. 8. — Radio Luminosity (1.4 GHz) vs. u-r color for a 
volume-limited subsample of the quiescent MAIN galaxy sam- 
ple. Since the fraction o f diluted AGN increases with redshift 
l)Revigl io & Hclfand 2002), the radio luminosity might be expected 
to increase if misclassified AGN dominate the radio emission. 

low most of the Ha emission. This is supported by some 
[Nil] emission even in the absence of Ha. The inset of 
Figure [SJd shows a zoomed- in view of this region for the 
high-SFR galaxies in the MAIN sample. Ha emission is 
just visible over the accompanying deep absorption fea- 
ture. Some may argue that the [Nil] /Ha ratio is more 
easily explained by a LINER (low ionization emission- 
line region) spectrum, and this may well be the case. 
However, there is also an obvious turn-up below ~ 3300 
A in the spectrum of the high-SFR MAIN galaxies, and 
no sign whatsoever of a similar trend in the low-SFR 
MAIN galaxies or LRGs. 

For both samples, the spectra with higher median 
SFRs show greater [Oil] emission a useful proxy fo r star- 
formation in the absence of Ha (|Kennicuttl H992) . (We 
should note, however, that lYan et al.l (|2006[ ) found [Oil] 
emission in red and post-starburst galaxies can be as 
strong as in star- forming galaxies.) We have attempted 
to quantify this result for the MAIN sample galaxies in 
Figure [7l The [Oil] Luminosities on the vertical axis 
were calculated by determining the equivalent widths and 
using the conversion to luminosity given in (Ke nnicutt] 
1992). They are only rough estimates, however, as we 
used magnitude in the (/-band rather than the B-band, 
and we used the median magnitude for the conversion of 
the composite equivalent width. For comparison, we also 
plot the results for the Star-forming sample. The [Oil] 
Luminosities for the MAIN sample do increase with SFR, 
but only weakly relative to the Star-forming sample. If 
we assume the [Oil] in the quiescent MAIN sample traces 
SF just as reliably as in the SF sample, and if we use the 
SF galaxies to calibrate an [OII]-SFR relation, the im- 
plied SFRs are all < 1 M yr -1 . There also appears to 
be a kink around 1 M Q yr , the same point at which 
the radio luminosities diverge, prompting speculation of 
a change in population at that point. Regardless, the 
weakness of the [Oil] dependence on D4000 SFR may in- 
dicate that both SF and AGN activity are contributing 
to the line emission, with the AGN component dominat- 
ing. Of course, why these galaxies have such weak line 
emission has been the question all along, and it pertains 
to both AGN and SF. 

The real question here is which component (AGN 
or SF) dominates the radio emission, and we do 



Fig. 9. — Radio Luminosity (1.4 GHz) vs. z-band Absolute 
magnitude for quiescent and star-forming galaxies. Arrows indicate 
upper limits. The star-forming galaxies have been split into 4 bins 
based on their fiber-corrected star-formation rates, and the median 
SFRs and [SFR ranges] are shown in the legend. We fit a line to the 
quiescent galaxy data; the least-squares fit gives log(L^) = (-0.92 
± 0.02)M Z + (7.5 ± 0.4). Additional error bars are overplotted on 
the star-forming galaxy points. These larger error bars correspond 
to 10% / 90% of the SFR dispersion for the galaxies within each 
bin. 

not believe that these composite optical spectra alone 

can answer that question. I n a recent paper, 

iMagliocchetti. Andreani. fc Zwaanl (|2007h use a sample 
of optically obscured, 24-/im-selected galaxies observed 
by Spitzcr in the First Look Survey (FLS) to con- 
clude that the 1.4 GHz luminosity is dominated by star- 
formation activity despite the presence of an AGN com- 
ponent. IRAC photometry sets the majority (66%) 
in a redshift interval of z = 1-3. The same con- 
clusion was drawn in Paper II for optically classified 
AGN/Composites in the local universe with SFR > 1 
M Q yr" 1 . 

If AGN activity is the source of the radio emission 
from the MAIN galaxy sample, then one possible ex- 
planation for the weak line emission is dilution. It has 
been shown that stellar continuum light of a host galaxy 
may dilute the n uclear emission lines of an AGN and hin - 
der its detection (iMoran. Filippenko. fc Chornockll2002f) . 



iReviglio fc HelfandT i 20071 ) present evidence that this ef- 
fect leads to an increasing number of misclassified weak- 
lined AGN with increasing redshift in the SDSS and other 
redshift surveys. We test for this effect in Figure[Hl where 
we take a volume-limited subsample in a narrow optical 
luminosity range and assess the effect of redshift on radio 
luminosity versus u-r color. Note that the magnitude cut 
of -22 < M z < -23 has preferentially selected the more 
luminous galaxies (in the radio as well as the optical). 
If increasing redshift bins show larger luminosities, we 
might conclude that AGN are contributing significantly 
to the radio emission. While there seems to be a slight 
trend, the effect is quite small and the error bars are 
large. The null result is therefore inconclusive. 

If the radio emission of the MAIN galaxy sample stems 
from star- formation, dilution might still be the culprit. 
To test this theory, we plot radio luminosity versus z- 
band absolute magnitude for the quiescent and star- 
forming galaxy samples in Figure[9] The SFR axis on the 
right-hand side of this figure (and those that follow) uti- 
lizes the radio-SFR calibration equation obtained in the 
previous section using the star-forming galaxy sample. 
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Fig. 10. — Radio Luminosity (1.4 GHz) vs. z-band absolute 
mag nitude for quiescent MAIN sample galaxies and LRGs. SFRs 
from Brinchmann ct al. ( 2004) arc ovcrplotted for each bin of both 
samples. 
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Fig. 11. — Same as Figure [9] but with UV-derived SFRs from 
GALEX ovcrplotted for each point as open squares. Though the 
UV-derived SFRs and the D4000 SFRs are discrepant by 1-2 orders 
of magnitude for the MAIN galaxies, we find them to be tightly 
correlated, arguing for a common origin. 

If we suppose for the moment that the radio emission 
of the Main sample quiescent galaxies originates with 
star-formation, then we can apply it to these galaxies 
as well (it obviously does not apply to the LRGs given 
our above results). We derived the mass (upper) axis 
from a linear least-squares fit to the plot of log(stellar 
mass) versus M z , where the masses came fr om the Stel- 
lar Mass Catalog (|Kauffmann et al.l [2003bh . The star- 
forming galaxy sample was first sorted by SFR and split 
into 4 subgroups to aid in the interpretation of the qui- 
escent galaxy SFRs. The median SFRs and ranges in 
median SFR for the bins in a particular subgroup are 
shown in the legend. The stacks were reduced to 1500 
objects per bin to allow for an adequate number of data 
points in each subgroup, and the resultant subgroups are 
plotted as different symbols for easy differentiation. The 
data for the LRG and MAIN samples were binned and 
stacked following our standard procedure, with the error 
bars representing ±1ct in median radio luminosity. We 
find that the quiescent MAIN sample follow a power-law 
fit of log(L^) = (-0.92 ± 0.02)M Z + (7.6 ± 0.5). 

While the Ha-based method does yield a large overall 
range of star-formation rates, the SF sample display a 
narrow range of SFRs for a given absolute magnitude of 



the host gal a xy. T his result was previously confirmed by 
iSalim et alJ £2007). If we assume no AGN contamina- 
tion for the moment, then using stacked radio luminos- 
ity as an indicator yields SFRs for the quiescent MAIN 
galaxies that are consistently 3-5 times lower than the SF 
galaxies. Keep in mind that we are not suggesting that 
the quiescent MAIN sample galaxies are intrinsically the 
same as the line-emitting SF galaxies, which we already 
demonstrated to populate a tight sequence in SFR-M Z 
space. But if we are to believe the trend demonstrated in 
FigurelH particularly for the low-L^j end, then one inter- 
pretation is that the strength of Ha needed for detection 
(and for subsequent SFR-estimation) is proportional to 
absolute magnitude. This interpretation implies a sort 
of Ha-dilution. To illustrate this point further, we have 
taken the star- forming galaxies in Figure [9] and ovcrplot- 
ted larger error bars representing the 10% / 90% range 
in SFR dispersion for the galaxies comprising each data 
point. 

Figure [10] uses the same axes to show how the radio- 
derived star-formation rates for th e quiescent galax- 
ies co mpare to those available from iBri nchman n et al.l 
in a different way than Figure [U The 8-point 
stars still represent the quiescent galaxies targeted in the 
MAIN Galaxy Sample, while the objects in the LRG sam- 
ple are shown as open circles in the upper left-hand re- 
gion of the plot. We have then overplotted (as crosses and 
squares, respectively) the median optically-derived SFR 
for each bin using the SFR axis on the right-hand side 
of the plot. For the quiescent LRG sample, the D4000 
SFRs are roughly two orders of magnitude lower than 
those inferred from the radio data; i.e., there is excess 
radio emission over that which is attributable to star- 
formation. This is again consistent with the hypothesis 
that the radio emission from these LRGs is dominated 
by AGN activity. 

The data for the MAIN sample quiescent galaxies 
show somewhat better agreement. The D4000 SFRs 
for the MAIN sample overlap the radio-derived SFRs 
at a few M Q yr" 1 and an absolute magnitude of M z 
= -22.5, though the two estimates show different slopes 
with absolute magnitude. If the quiescent MAIN sam- 
ple are star-forming, this discrepancy could simply be 
explained by the different methods used to calculate 
SFR between SF and "Unclassifiable" galaxies. If the 
emission is dominated by an AGN contribution, alter- 
natively, it must scale with galaxy stellar mass to ex- 
plain the observed correlations. This is not surprising 
for star formation in a homogenous sample of galaxies, 
but one might argue that this is expected for AGN as 
well, given the well-known scalin g relation between black 
hole mass and bulge luminosity ( Kormendv fc Richstonel 
U995tlFerrarese fc M crritt 2000). Our own results for the 
quiescent LRGs may even support this claim. 

As we have discussed above, radio emission is an in- 
direct indicator of star-formation. The UV continuum, 
on the other hand, is thought to be directly p roportional 
to th e photospheric emission of young stars (jKennicuttl 
1998). Still, it is not free from problems of its own, most 
notably significant dust obscuration. Nevertheless, as 
another c omparison, we ove rplot UV-derived SFRs from 
GALEX (jSalim et alJ l2007f ) in Figure Q3J corrected to 
our IMF. As our samples do not completely overlap, we 
have simply taken the median of the UV SFRs available 
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for each bin. For our MAIN sample quiescent galaxies, 
there were generally ~ 300 galaxies (per 3000-galaxy bin) 
with UV SFRs, and we find these galaxies to have fairly 
representative radio luminosities (i.e., there is no selec- 
tion bias between our samples). None of the quiescent 
galaxies targeted as LRGs had UV-SFRs available from 
this catalog. 

The se UV-based SFRs were obtained by ISalim et"aT1 
(2007) via SED-fitting of emission f rom the UV to the 
z-band with various models from the iBruzual fc Chariot! 
(2003) population synthesis code. Because the UV lumi- 
nosity imposes the tightest constraint on the SFR, the 
authors refer to these SFRs as "UV-derived" . Their 
quoted rates are averaged over a span of 100 Myr 
and constitute the first test of the two-component 
dust a ttenuation model proposed by I Chariot &; Fall 
(2000]). ISalim et alJ (l2007h co mpare these rates w i th the 
optically-derived rates for the iBrinchmann et all (|2004f ) 
star- forming sample and report good overall agreement, 
with some systematic deviation due to differences in 
modeling dust attenuation. It is therefore not surpris- 
ing that we also find good agreement in the star-forming 
sample, as they are derived from the same catalog and 
served as our radio-SFR calibrators. 

The interpretation of the UV and radio-derived SFRs 
for the MAIN sample quiescent galaxies is, once again, 
less straightforward. The UV-derived SFRs are be- 
tween 2-3 orders of magnitude lower than the SFRs de- 
rived from radio-stacking. Given the rough agreement of 
the radio and D40 00 SFRs, it is then no surprise that 
ISalim et alJ (|2007f ) find a ~ 2 order-of-magnitude dis- 
crepancy for the UV versus optically-derived SFRs of 
their "No Ha" sample, which comprise roughly ~ 10% 
of our total MAIN quiescent sample. They conclude that, 
since these galaxies fall largely on the red sequence, the 
UV SFRs are more realisti c. Another of their argum ents 
is that such high SFRs as IBrinchmann et ahl (|2004h es- 
timate would result in a large fraction of their "No Ha" 
galaxies having detectable Ha. On the contrary, we have 
shown in Figure[9]that this result could also be explained 
if the strength of Ha required to constitute a "detection" 
is simply a function of optical luminosity, with the levels 
predicted for these quiescent galaxies falling below the 
detection limit for a given absolute magnitude. 

The strongest argument for the UV-SFRs is that the 2- 
3 orders of magnitude discrepancy between the radio and 
UV would require an underestimation of dust extinction 
of 5-8 magnitudes. This is, admittedly, quite large, but 
perhaps the disagreement between SFR estimates could 
be explained by the location of the star-formation; If the 
star-formation is dominated by nuclear star-formation, 
it might be easier to hide both the line-emission and the 
UV flux from observation. This is speculative, but we 
would like to emphasize again that the radio emission 
and D4000 index are two independent checks on the pres- 
ence of star-formation, and they are consistent with one 
another to within a factor of ~ 3. While they both can 
be higher than the true SFR for completely different rea- 
sons, we feel this is a coincidence that requires further 
explanation before the possibility of SF in these galaxies 
can be thrown out altogether. We will, however, freely 
acknowledge that if the GALEX SFRs are right, then 
for sure the radio-derived SFRs are due to AGN activity. 
Overall, whether this discrepancy is ultimately the result 



of an AGN component to the radio emission or a poor 
dust-attenuation model is not a question we feel can be 
definitively answered here. 

5. CONCLUSIONS 

We have investigated the radio emission from a sam- 
ple of ~ 185,000 optically-selected quiescent galaxies 
from the Sloan Digital Sky Survey. By stacking FIRST 
cutouts centered on the galaxy positions, we reach me- 
dian radio flux densities down to the 10s of /LtJy. We 
calibrate the relation between 1.4 GHz radio emission 
and optically-derived SFR using a sample of star-forming 
galaxies, allowing us to look for star-formation in their 
quiescent counterparts. The main findings of this study 
are summarized below. 

• We find a tight correlation between radio luminos- 
ity and SFR for the star-forming galaxy sample, 
with Lx.iGHz going as (SFR) 137 . 

• The quiescent galaxy sample consists of two sub- 
groups: those targeted as part of the MAIN Galaxy 
Sample, and those targeted as part of the Luminous 
Red Galaxy (LRG) Sample. The former show par- 
tial overlap with the star-forming sample on a plot 
of \jr vs. optically-derived SFR, while the latter 
show a relatively constant, high value of median 
radio luminosity that does not vary with SFR. 

• For the LRGs, the correspondingly high rates of 
radio-derive d star-formation conflict w ith those es- 
timated by IBrinchmann et aTl (|2004h via calibra- 
tion of the 4000-A break by two orders of magni- 
tude on a plot of derived-SFR vs. M z . This re- 
sult, combined with the relative constancy of L/j 
with D4000 SFR, suggests that the radio emis- 
sion of these quiescent LRGs derives not from star- 
formation but from an AGN component. This im- 
plies the possible presence of radio jets and a mech- 
anism for low-level feedback in these galaxies. 

• The MAIN sample galaxies, on the other hand, 
exhibit a L^-SFR trend which intersects the SF 
sample around 1 M Q yr _1 and deviates by no 
more than a factor of 3 over the whole range of 
SFRs. This agreement as well as the dependence 
of Ln with SFR suggests that the L^ might arise in 
part from star-formation (particularly in so far as 
the AGN of Paper II show no such dependence.) 
If we assume the radio emission of these galax- 
ies originates from star-formation, the deviation 
at high/low values of SFR might then be due to 
the different diagnostic measures used to estimate 
SFR for the two samples (i.e., 4000-A break vs. 
emission-lines). However, we cannot rule out the 
possibility that the radio emission may have a large 
or even dominant AGN component, with a Lr- 
mass relation causing the observed correlation. 

• Composite spectra of the MAIN sample show 
Balmer absorption strong enough to plausibly swal- 
low any Ha emission present. This is supported 
by some [Nil] emission even in the absence of Ha. 
The [Oil] emission, however, increases only slightly 
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with D4000 SFR, giving reason to doubt the valid- 
ity of D4000 (and thus radio emission) as a proxy 
of star-formation in these galaxies until we under- 
stand why the line emission is so highly-suppressed. 

• The MAIN quiescent sample shows a correlation 
between radio luminosity and z-band absolute mag- 
nitude of the form: log(Li. 4 ) = (-0.92 ± 0.02)M Z 
+ (7.6 ± 0.4). The SFRs implied by radio emis- 
sion reach above 10 M Q yr _1 , yet lie just below the 
objects classified as "star-forming" for a given ab- 
solute magnitude. If we assume the radio emission 
stems from SF, one possible interpretation is that 
the strength of the optical emission lines needed to 
constitute a "detection" varies with galaxy lumi- 
nosity, and using our radio-stacking method allows 
us to probe low levels of star-formation despite SF 
dilution in the optical. If we assume the radio emis- 
sion stems from AGN activity, then both the D4000 
and radio-derived SFRs are too high (for separate 
reasons), and the observed trend implies a relation 
between AGN radio luminosity and galaxy mass. 

• UV-derived SFRs from GALEX are 2-3 orders of 
magnitude lower than the SFRs predicted by radio- 
stacking for the MAIN quiescent galaxy sample. 
This indicates possible AGN contamination in the 
radio emission of the quiescent galaxies, incom- 
plete modeling of the dust attenuation for the UV 
emission, or both. The different processes at play 
obviously need to be better understood before we 
can choose the most useful SFR diagnostic tool for 
these optically-quiescent galaxies. 

In short, we use a novel stacking technique to provide 
strong evidence that quiescent LRGs harbor AGN. This 
result has clear implications for galaxy formation sce- 
narios, providing observational support to the existence 
of a mechanism for low-level feedback in these galaxies. 
The results for the quiescent MAIN sample galaxies are 
more nebulous, with the discrepancies between the dif- 
ferent SF indicators clearly highlighting how much we 
still have to learn about these persistently unclassifi- 
able galaxies. Overall, we believe the technique utilized 
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